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Abstract 

We construct the sphaleron solution, at zero and finite temperature, in the 
Minimal Supersymmetric Standard Model as a function of the supersym- 
metric parameters, including the leading one-loop corrections to the effec- 
tive potential in the presence of the sphaleron. At zero temperature we have 
included the one-loop radiative corrections, dominated by the top/stop sec- 
tor. The sphaleron energy -Emssm mainly depends on an effective Higgs 
mass = lim m/^, where m/j is the lightest CP-even Higgs mass and 

rriA the pseudoscalar mass. We have compared it with the Standard Model 
result, with = mf^ , and found small differences (1-2%) in all cases. 
At finite temperature we have included the one-loop effective potential im- 
proved by daisy diagram resummation. The sphaleron energy at the critical 
temperature can be encoded in the temperature dependence of the vacuum 
expectation value of the Higgs field with an error < 10%. The light stop 
scenario has been re-examined and the existence of a window where baryon 
asymmetry is not erased after the phase transition, confirmed. Although 
large (low) values of rrih {rriA) are disfavoured by the strength of the phase 
transition, that window (along with LEP results) allows for rrih < 80 GeV, 
niA > 110 GeV and At < 0.4 tuq. 
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1 Introduction 



't Hooft observation [0] that baryon, and lepton, number is not conserved at the quan- 
tum level in the Standard Model (SM) of electroweak interactions, due to the presence 
of an SU(2) triangular anomaly, opened up the exciting possibility of baryon number 
non-conservation at low-energy. Even if the probability for AB ^ processes at zero 
temperature is exponentially suppressed by a factor ~ exp(— 167r^/(7^), where g is the 
SU(2) gauge coupling, such suppression factor is absent at high temperatures 
and, more precisely, at temperatures of the order of 100 GeV, triggering the hope of 
generating the baryon asymmetry of the Universe |§] at the electroweak phase transi- 
tion 0. 

The existence of sphalerons, static and unstable solutions to classical field equations 
in the SU(2) gauge theory with Higgs fields in the fundamental representation, has been 
known since long ago [|1^. A key observation was that of Klinkhamer and Manton 



who interpreted the sphaleron solution as a saddle point of the energy barrier separating 
gauge-inequivalent classical vacua. Detailed numerical calculations of the sphaleron 
solutions and energy in the SM were subsequently examined by many authors 0, ^, ^. 
The precise value of the sphaleron energy is of the utmost importance for mechanisms 



aiming to explain the generation of the baryon asymmetry of the Universe ||ll||-|ll4|. 
In particular the sphaleron transition rate for temperatures below the phase transition 
temperature 

r,ph~exp<^ — - — ^ (1.1) 
along with the bound obtained by comparing Eq. (|L1|) with the expansion rate of the 



Universe at the temperature of the electroweak phase transition [0, [15 



> 45 (1.2) 



where Tb denotes the transition, or bubble formation temperature, allows to put upper 
bounds on the Higgs boson mass and thus discard models of electroweak symmetry 



breaking on the basis of baryon asymmetry generation. 

Even though the SM has all the required properties for the generation of the baryon 
asymmetry (CP violation, baryon number violating processes and non-equilibrium pro- 
cesses generated by a first-order phase transition) they are not strong enough to gen- 
erate the necessary amount of baryon asymmetry. In particular, the CP-violating 
Cabibbo-Kobayashi-Maskawa phases strongly restrict the possible baryon number gen- 
eration JlBl, though the latter remains as a very controversial subject. Even disregarding 



this point, the phase transition is not strongly first order enough and, moreover, con- 
dition (|1.2| ) translates into an upper bound on the Higgs boson mass [|l^, [l^, |19| which 
is well below the present experimental lower bound of 65 GeV. 

Imposing the requirement of baryon asymmetry generation at the electroweak phase 
transition translates into the requirement of new physics at the weak scale. As a 
matter of fact, when extending the SM at the weak scale and probing the baryogenesis 
capability of the extended model the detailed calculation of the sphaleron energy should 
be an essential piece of information. It has been proved that there are extra sources of 
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CP violation, suitable for baryogenesis, in minimal extensions of the SM, in particular in 
the SM with a gauge singlet [^0|, |2ll] (with CP-violating non-renormalizable couplings 
arising from integrated out extra fields), and in two-Higgs doublet models ||22| (with 
extra phases in the Higgs sector), and that the phase transition is strongly first order 
enough not to wipe out, after the phase transition, the previously generated baryon 
asymmetry |]23| , ^ for values of the Higgs mass beyond the experimental bounds. In 
both cases the sphaleron energy has been computed, Refs. |^ and |^ respectively, and 



can support the baryogenesis achievements of the corresponding models. However, the 
sphaleron solutions have been evaluated in these models using the tree-level potential 
at zero temperature, so that the sphaleron energy at the phase transition temperature 
relies on the assumption of the scaling law 

i?sph(T)=i?sph(0)^m (1.3) 



which has only been proved to hold approximately in the SM case , a similar calcu- 
lation being absent in the above extensions of the SM. 

The possibility that new physics at the electroweak scale be provided by the Minimal 
Supersymmetric Standard Model (MSSM) is a very appealing one on the theoretical 
front, since it can technically solve the hierarchy problem. The MSSM is also the natural 
candidate for an effective theory from a more fundamental theory valid at the high scale 
(^ Mpi), as e.g. string theory. Finally LEP electroweak precision measurements predict 
in the MSSM the unification of gauge couplings at a scale ~ 10^^ GeV, which is the 
only present 'evidence' for Grand Unification and new physics beyond the SM. All 
of that has generated a plethora of experimental searches for new physics within the 



MSSM framework in present and planned colliders |2^, 0. Exploring the capability 
of the MSSM of generating the observed baryon asymmetry of the Universe at the 
electroweak phase transition is therefore of the highest interest, mainly in order to 
confront the experimental searches of supersymmetry with possible regions of the space 
of supersymmetric parameters where baryogenesis might be possibly produced. In this 
sense new sources of CP violation are present in the MSSM [0, which can serve 
to overcome the strong SM suppression factors mentioned above []TB|. On the other 
hand, the strength of the first order phase transition has been extensively studied in 
the MSSM [^-|^, and proved that the weak phase transition in the SM (driven by 
the gauge coupling) can be strengthened in the MSSM (driven by the top Yukawa 
coupling) in the presence of light supersymmetric partners of the top quark (stops) and 
small values of tan /5 . However a detailed study of the sphaleron solutions and the 
sphaleron energy in the MSSM was missing. 

In this paper we will construct the sphaleron solutions in the MSSM at zero and 
finite temperature, including in all cases the leading one-loop radiative corrections. The 
planning of this paper is as follows: In Section 2 we will construct the sphaleron so- 
lutions and the sphaleron energy in the MSSM at zero temperature, including the full 
set of one-loop radiative corrections to the effective potential in the presence of the 
sphaleron Higgs. They are dominated by the top Yukawa coupling, and we will consis- 
tently neglect the bottom Yukawa coupling, a reasonable approximation for tan (3 ^ 15, 
which are the values we will consider. In fact, as we will show, only for tan/3 < 3 
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the baryogenesis scenario is feasible. We have compared in every case the sphaleron 
energy in the MSSM, with arbitrary sets of supersymmetric parameters predicting the 
supersymmetric and Higgs mass spectra, with the sphaleron energy in the SM with an 
effective Higgs mass. The deviations are typically ^ 1.5%. In Section 3 we will con- 
struct the sphaleron solutions and energy in the MSSM at finite temperature. We have 
included the one-loop effective potential at finite temperature, in the presence of the 
sphaleron Higgs, improved by the resummation of daisy diagrams. We have checked 
that the approximation given by Eq. (|1.3| ) is accurate with an error less than a few per- 
cent for the cases where the phase transition is weak. For the cases studied in Section 
4, where we have applied the results of Sections 2 and 3 to the recently proposed case 
of light stops, the error can be as large as ~ 10%. We have qualitatively confirmed the 
results of Ref. |^6|]. In particular we have found a baryogenesis window for rrij < m^. 



nih 80 GeV, tan/5 < 3 and uia ~ HO GeV, although not all the previous inequalities 
can be simultaneously saturated. Finally Section 5 contains our conclusions. 

We have adopted, throughout this paper, the approximation g' = 0, where g' is the 
U{1)y gauge coupling, in order to use a spherically symmetric ansatz for the sphaleron. 
The corrections of 0{g''^) to the sphaleron energy have been evaluated in Ref. 0, com- 
puted in Ref. for the SM sphaleron, and proved to be negative and < 1%. Also 
we have kept in the effective potential only one-loop corrections in the presence of the 
sphaleron Higgs, which are dominated by the top/stop sector and the top Yukawa cou- 
pling. We have neglected one-loop radiative corrections in the presence of the sphaleron- 
W . This procedure should provide a rather accurate one-loop approximation because 
the leading finite temperature 0{T'^) corrections to the effective potential in the pres- 
ence of the sphaleron-iy vanish [} Subleading finite temperature corrections to the 
effective potential are suppressed by powers of n > 2, and can therefore be ne- 
glected as compared to those proportional to powers of the top Yukawa coupling. The 
approximation of considering only radiative corrections to the effective potential (or 
effective action) in the presence of the sphaleron Higgs, and neglecting those generated 
in the presence of the sphaleron-iy, is consistent only in cases where the former are 
dominated by a large coupling (e.g. in our case by the top Yukawa coupling) while the 
latter are always provided by extra powers of the gauge coupling constants. In other 
words, the approximation should work in theories where putting g = Q m radiative 
corrections is a good approximation. An example of this kind of theories is the MSSM, 
where radiative corrections are dominated by the top Yukawa coupling. Unlike the 
MSSM, the finite temperature radiative corrections in the SM are dominated by the 
gauge coupling and making this approximation would miss some terms similar to those 
considered [p^ . 



^This is due to the fact that the sphaleron solution is non-zero only for transverse degrees of freedom 
whose (magnetic) Debye mass vanishes to all orders in perturbation theory This is not the case for 
longitudinal degrees of freedom, so that a non-vanishing component along Wt^ would trigger, through 
the one-loop electric mass, sizeable corrections. 
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2 Sphalerons at zero temperature 



In this section we will compute the static unstable solutions of classical equations of 
motion in the MSSM. As we have said in Section 1 we will work in the approximation 
of taking g' — so that the U{1)y gauge field can be consistently set to zero. This 
will allow a spherically symmetric ansatz, as follows. The lagrangian density for the 
SU (2) gauge fields and the Higgs system 



Ho 



Hi 
Ht 

^2° 



(2.1) 



is given by 



>C = --W^.W''^' + {D^H,)"^ {D'^H,) + {D^H^y {D^^H^} - V,s{Hu H^) 
where the SU{2) field strength is defined as, 

the covariant derivatives are 



2 



(2.2) 



(2.3) 



(2.4) 



cr" being the Pauli matrices, and Ves is the effective potential where we will consider 
one-loop corrections. 

We can expand the effective potential (at zero temperature) as 



V,s = Vo{H,,H2)+V^{H,,H2) + 



(2.5) 



where Vq is the tree-level potential and Vi contains the one-loop radiative corrections. 
They are determined from the supersymmetric structure of the MSSM, with superpo- 
tential 

W = htQL ■ H2UI + At//i • H2 (2.6) 
and from the soft-breaking terms. In particular, the tree-level potential is given by 



Vq{H^, H2) = ml hIHi + ml HIH2 + ml{H^ ■ H2 + h.c.) 

(hIh, - H^H^y + 4 (HIH2) (hIh,) 



(2.7) 



and the one-loop corrections, in the 't Hooft-Landau gauge and in the DR renormal- 
ization scheme, by 



V^{H^,H2)=j:-^2^t(H„H2) 



log 



mj{H^,H2) 3 
g2 2 



(2.8) 
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where Q is the renormahzation scale, that we are taking, for definiteness, as = m: 
m]{Hi, H2) is the field dependent mass of the i^^ particle in the background Hi, H2, 
and rii is the corresponding number of degrees of freedom, which is taken negative for 
fermions. The relevant degrees of freedom for our calculation are the gauge bosons 
(W, Z), the top quark (t) and its supersymmetric partners (ti, ^2), with 



n 



t2 



Tit = -12, nw = 6, nz = 3, rij^ 
Choosing the temporal, radial gauge 

Wo = 0, x,W,^ = (i = 1 - 3, a = 1 - 3) 
we can write the static, spherically symmetric ansatz 

2/(r 



6. 



(2.9) 



(2.10) 



W7(x) 



Hi{x) = hi{r) i 
H2{x) = h2{r) i 



gr- 

a ■ X 

r 

a ■ X 



1 




1 



(2.11) 



for the sphaleron in the MSSM, where = + + z^. 

There, however, remain time-independent gauge transformations which can be used 
to rotate the background fields in ( |2.11 ) into a more convenient form for our calculation. 
In particular, the gauge transformation 0] 



H (J ■ X 



U (x) = exp 



transforms the ansatz ( p.llj ) into the background fields 



(2.12) 



Hi{x) = hi{r) 
H2{x) = h2{r) 



1 



1 



(2.13) 



which should smoothly approach a vacuum in the unitary gauge as r ^ 00. 
In the presence of the background ( |2.13| ) the tree-level potential reads as 

Vo{hi, h2) = ml hl{r) + mj hl{r) + 2 mlhi{r)h2ir) + ^ \hl{r) - hl{r)'^ 



(2.14) 



while the one-loop corrections are given by Eq. ( ^.81 ) with background dependent masses. 
The relevant masses are the top quark mass, given by 



miih) = hthiir) 



(2.15) 
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the gauge boson masses, by 

,2 



and the stop squared mass matrix given by 



with entries 



2 2 

m~ m~ 

2 2 

m~ m~ 



(/i) = m| + m2(/i) + {g^ [h\{r) - hl{r)] 



(2.16) 



(2.17) 



m~ (/i) = + mf{h) 



(2.18) 



m^- ^(/i) = /it (A/i2(r) + /i/ii(r)) 



the mass eigenstates being defined by the diagonahzation of matrix ( [2. IT] ), 



m~ (h) + 171% (h) 
'fn-r \h) = ± 



m~ ih) — m~ ih) 



- 2 




+ 


m~ [K) 

iL,R. 



(2.19) 



with h~ {hi, /i2). 

By minimizing the effective potential ( |2.5|) with respect to hi,h2, and imposing the 



minimum of the potential at (fi,f2), witht> = yt>f + f I = 171.1 GeV, and tan/? = V2/V1 
fixed, we can eliminate ml and 1712 in favour of the other parameters of the theory, 
as 



35, 41 



-mgtan/?- ^cos2/3-^ 



647r2 



-m 



cot/5+^cos2/3-y At 



h={vi,V2) 



h={vi,V2) 



(2.20) 



while ml can be traded in favour of the one-loop corrected squared mass m\ of the 



CP-odd neutral Higgs boson, as ^ 
777.3 ~ ~''^A sill P cos /3 



\og{m\/Q^)-l 



t2 



log(m|/Q2) _ 1 



32n^ml, sin^ /? 



2 2 
m ~ — m- 



t2 



(2.21) 



/l=(l'l,D2) 



The energy functional corresponding to the ansatz of Eq. ( |2.13| ) can be written as: 



static 



47r I drl^ 



(drhf + -hiii-fy 



{drh2f + -hl{l-ff 



(2.22) 
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where we have subtracted the vacuum energy. The sphaleron is described by the func- 
tions f{r), hi{r) and h2{r) which are the solutions to the equations of motion 



(9,2 + \) h^ir)-^hir)il-f{r)y 



h2ir)--h2{r){l-f{r)y 



ldV,s{hi,h2) 
2 dhi 

2 dho 



dl fir) - 4/(0(1 - /(r))(l - 2/(r)) = -^(1 - fir)) [hUr) + hl{r) 
subject to the boundary conditions, 

/ii(0) = /i2(0) = /(0) = 

/ll(oo) = fi, /l2(oo) = V2, /(OO) = 1 



(2.23) 

(2.24) 
(2.25) 



(2.26) 



which guarantee that the solutions of ( p.23[ - pl25D have finite energy and the correct 
vacuum behaviour when r — 00. 

In this section we have introduced the set of supersymmetric parameters which the 
sphaleron solutions and sphaleron energy depend upon. They are tan/5 = V2/V1, m^, 
the mass of the CP-odd Higgs boson, fi the supersymmetric higgsino mass, mq and 
mfj, the soft supersymmetry breaking squared mass terms for the left-handed doublet 
and right-handed singlet squarks 



Qi 




(2.27) 



and At the soft breaking term for the third generation of squarks corresponding to the 
tri-linear coupling in the superpotential (|2.6| ). We will consider tan (3 < 15 and therefore 
we can neglect the bottom Yukawa coupling hh, and thus the corresponding tri-linear 
soft term Ah will have no influence on the result, and can be neglected too. As for 
the top Yukawa coupling, we will fix rrit = 175 GeV throughout this paper. This 



value of rrit is inside the experimental range [^2|, ^ and we have decided to fix it in 
view of the large number of parameters in the model. Variations of the results in this 
paper corresponding to variations of rrit inside the experimental error band will not 
dramatically modify the final conclusions. 

In Fig. |1] we have illustrated the functions hi{r)/v, h2{r)/v and /(r) [solid lines] 
for a particular value of supersymmetric parameters: tan/5 = 1.5, mq = 500 GeV, 
rriji = 100 GeV and m^/ = At = fi = 0. In Fig. |^ we plot the sphaleron energy given 
by Eq. (|2.22|) [in units of m^/aw, oiw = 5'^/4vr] as a function of tan/? for different 
values of the supersymmetric parameters, where we have varied mu = — 400 GeV and 
rriA = 100 — 500 GeV. Comparison of curves (a) and (b), and (c) and (d), shows that 
the variation of i^sph with m^, within the considered range, is tiny, while comparison 

^Notice that mt is the MS runmng top mass determined at the scale Q ~ Mt. It is related to 
the pole top mass Mt by the QCD correcting factor mt = Mt/ [1 + 4as(Mt)/37r]. The on-shell mass 
mt = 175 GeV corresponds to Mt = 183 GeV. 
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of curves (a) and (c), and (b) and (d), exhibits the influence of the parameter mjj in 
the value of -Egph, being < 3%. We can see from Fig. ^ that i^sph varies ~ 10% over the 
whole range of tan /3. 

A similar exercise to that performed in Fig. ^ has been done in Fig. |^ exchanging 
the variable tan/3 with the variable m^. In Fig. ^ we plot i^sph as a function of tua for 
different values of the other supersymmetric parameters where we have varied tan /3 = 
2 — 15 and mjj = — 400 GeV. The large variation with tan/? is explicit by comparing 
curves (a) and (b), and (c) and (d). The variation with uiu, comparison of curves 
(a) and (c), and (b) and (d), being smaller. Finally the variation with is tiny, as 
previously noticed, for the considered range of supersymmetric parameters. 

The behaviour of i?sph in Figs. |^ and ^ as a function of the different supersymmetric 
parameters can be understood in a simple way. The Higgs boson mass spectrum is 
determined in the MSSM from the tree-level potential (|2.7| ) and the radiative corrections 
In particular the CP-even Higgs boson masses rrih^H can be determined from the 



supersymmetric parameters [44]- 



mh,H = mh^nimA, tan/?, ttiq, niu, At, /i, . . .) (2.28) 

where the first two parameters determine the tree-level expression while the rest enter 
into the radiative corrections. Now, the field in the {hi,h2) plane which acquires a 
vacuum expectation value (VEV) is along the direction hi cos [3 + h2 sin (3, and the 
squared mass (curvature) of the potential along that direction is 



i->i I ' ^ 

m 



= sin2(a -/3)ml + cos^{a - (3) mjj (2.29) 

where a is the mixing angle in the Higgs sector, whose determination, as for the masses 
rrih and mn in (|2.28 ), involves radiative corrections. The simplest way of evaluating 



( p. 291 ), with radiative corrections provided by ( |2.8| ), is by observing that it does not 
depend on m^, and therefore, using the limit (when oo), sin^(a — /?)—*> 1, one 

obtains 

mf = lim ruh . (2.30) 

m^i— >oo 

Since the sphaleron solution asymptotically points towards hi cos f3+h2 sin f3 one expects 
that Esph will mainly depend on g. Notice that, from (^^301) , mj^ and nih will 
coincide only for very large values of m^. This behaviour is exhibited in Fig. ^ where 
we plot i?sph as a function of m'f for the same cases as in Fig. ||, and compare it with 
Ef^ for mf^ = m'f . We can see that the spreading is small (~ 1%) and due to the 
fact that the direction which does not acquire a VEV also contributes to the sphaleron 
functional. 

Finally the dependence of i^sph with the stop mixing parameters in Eq. ( p.l8| ). At 



and /i, are shown in Figs. ^ and |^. We can see from Fig. |^ that the variation of E'sph 
with At is < 3% for generic values of the supersymmetric parameters, while the variation 
with /i is smaller, as can be seen from Fig. |^. This fact can be understood from the 
fact that the effective stop mixing parameter is provided by 

At = At + f^/tanp (2.31) 



^We thank J.R. Espinosa for making this observation to us. 
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so that the effect of the /i parameter is suppressed by 1/ tan/3. For large values of tan/3, 
Esph does not depend on the fi parameter, as can be seen from curve (b) in Fig. 



3 Sphalerons at finite temperature 

In order to introduce finite temperature effects we have to modify the effective potential 
of Eq. ( |2.5|) by adding the thermal corrections. We will consistently work at the level 
of the one-loop thermal corrections improved by the resummation of daisy diagrams. 
The effective potential can be written, to this level of approximation, as 

Kff = Vo{h) + Vi{h) + AV,{h,T) + AKlaisy(/i,T) (3.1) 

where Vq and Vi are the tree-level and one-loop potentials at zero temperature, given by 



Eqs. (|2.14|) and (|2^ ), respectively, AVi the one- loop correction at finite temperature, 
and AVdaisy the resummation of daisy diagrams. They are given by: 



AV,{h,T) = f^ 



2^2 



(3.2) 



where rii is given in Eq. (|2.9| ), the masses mi{h) are defined in Eqs. ( p. 151 - ^^191) , and the 
thermal function Jj = J+(J_), if the i^^ particle is a boson (fermion), with 



J±(y2) = dx log (^1 T e-v^^ 



As for the term AVdaisy, it is given by 

AKlaisy(/^,T) = --^ 



m\{h,T) - m\{h) 



B 



(3.3) 



(3.4) 



where the sum is extended to scalar bosons and the longitudinal degrees of freedom of 
gauge bosons, with ub given by Eq. ( |2.9|) for ti^2 and 



nwr 



2, nzr 



1. 



(3.5) 



The thermal masses m^(/i,T) are obtained from m\{h) by adding the leading T de- 
pendent self-energy contributions, which are proportional to T^. In particular, the stop 
squared thermal mass matrix is given by 



_2 / m\ + (T) 

Mj=( *^ 2 ^ 2 

\ tLR tu tR 



tLR 



(3.6) 



The self-energies are given by 



^1r(T) 



12 



'a! + 1^ 7. + sin^ /3 + cos' /3 0^ /i? + 



tR 



1 ^2 



^9s + U^7 + sin' p + cos' peA)h 



tL 



(3.7) 
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where gs is the strong gauge couphng. Only loops of gauge bosons, Higgs bosons and 
third generation quark/squark have been included, assuming that all the remaining 
supersymmetric particles are heavy and decouple by Boltzmann suppression factors. In 
particular the gluinos, if light, would provide contributions to the self-energies which 
are given by 2/9 ^^t2. They constitute the main missing contribution to the thermal 
masses and would weaken the strength of the phase transition. We are considering 
them heavy as in Ref. |3^. We also have introduced exphcit step-9 functions for the 
contribution of left-handed and right-handed stops and pseudoscalar bosons, with the 
convention: Bx = 1 [0] for mx ^ it T [nix ^ vr T] . 
The thermal Wl-t^slss is given by 

m^J/i,T) = ml,{h) + UwAT) (3.8) 

where the self-energy is 

U^,{T) = ^-g'T' (3.9) 

and of course, since we are taking g' = 0, rhzj^ih, T) = m^^^{h, T). 

Up to now we have considered one-loop corrections to the effective potential in 
the sphaleron background (ansatz) corresponding to diagrams with hi{r) and h2{r) as 
external legs. Strictly speaking we should also consider diagrams with W°' as external 
legs. They should contribute to the effective action as 



-\K^:^^t + --- = -\9^T^W^WS + --- . (3.10) 



While the leading term in Eq. ( p.lOj ) cancels for the sphaleron ansatz of Eq. ( p. 11 ), the 



ellipsis denotes terms which are suppressed by powers of the gauge coupling constant 
and by inverse powers of the temperature and we will not consider them explicitly in 
this paper. As previously stated this procedure is self-consistent since in the MSSM 
the radiative corrections are dominated by the top Yukawa coupling and the top/stop 
sector. 

Now, at finite temperature, the energy of the sphaleron, the equations of motion 
and the boundary conditions for functions hi{r,T), h2{r,T) and f{r,T), are given 
by Eqs. { ^72"^^ ), where Vesihi,h2,T) is provided by Eq. (^ and vi = vi{T), 
V2 = V2{T) is the minimum of the effective potential (|3.1|) at the temperature T. We 
define the critical temperature Tc, as in Ref. as the temperature at which the 
determinant of the second derivatives of Ves{h,T) at h = vanishes i.e. 



det 



d^V,s{h,T,) 



' J J h\=h2=0 



0. (3.11) 



Explicit formulae for computing the critical temperature can be found in Ref. 35 1. For 



illustrative purposes we have plotted in Fig. ^ the functions hi{r,Tc)/v, h2{r,Tc)/v and 

^The temperature Tc, as defined in ( |3.1l| ), is usually called lower metastability temperature or 
lower spinodial decomposition point. We thank M. Shaposhnikov for pointing out this to us. We call 
it critical temperature just for simplicity. 
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f{r,Tc) (dashed lines) for tan/? = 1.5, = 100 GeV, mq = 500 GeV, niu = and 
At = fi = 0. 

The T-dependence of i^sph is shown in Fig. |^ for three typical sets of supersymmetric 
parameters (solid lines). In all cases the pattern of curves is similar to that shown in 
Ref. for the SM case. In particular Esph decreases with increasing temperatures 



and sharply goes to its minimal value at the temperature where the local minimum we 
are considering disappears. We have also probed the approximated scaling law, 

KltiT) = i?sph(0)^ (3.12) 



where v{T) = ^Jvf{T) +V2{T), fi(T) and V2{T) being the vacuum expectation values 
of the fields hi and /12 at finite temperature. We have plotted in Fig. |^ E^^^(T) (dashed 
lines) for the three cases. We have found that Eq. ( |3.12| ) is an excellent approximation, 
for all values of supersymmetric parameters used in Fig. |^, with an error < 3%. In fact 
we have found that the scaling law ( |3.12| ) is not controlled by any particular combina- 
tion of supersymmetric parameters, but by the strength of the phase transition. The 
weaker the phase transition the better the scaling law. In fact for a second order phase 
transition the scaling law would be exact as happens in the SM. The reason of the little 
departure between solid and dashed lines in Fig. ^ is because of the weakness of the 
phase transition for the particular values of the chosen parameters. In cases where the 
phase transition is stronger, as we will discuss in Section 4, the departure corresponding 
to the approximation of Eq. ( p.l2| ) is greater. 

Finally we have plotted in Figs. | and ^, for the cases considered in Figs. ^ and |^, 
Esph(Tc)/Tc as a function of tan P and m^, respectively. We can see that the dependence 
of Esph{Tc) on the supersymmetric parameters, is qualitatively similar to that found for 
-E'sph(O) in Figs. |^ and ^. We can see from Figs. ^ and |^ that for the sample of generic 
cases considered up to now the phase transition does not satisfy condition (|1.2|) and 
therefore is not strong enough first order to generate baryon asymmetry. This problem 



was taken care in Ref. where the order of the phase transition was considered and a 
'non-generic' region of the space of supersymmetric parameters was proposed to trigger 
a strong first order phase transition. This case is characterized by light right-handed 
stops, heavy left-handed stops and low values of tan f3 and will be analysed in the next 
section. 



4 The case of light right-handed stops 

In the previous sections we have studied the sphaleron solutions and energy in the 
MSSM for generic values of supersymmetric parameters, as was done in analyses of the 
phase transition in the MSSM in the previous literature P3-P3[. Our results on the 



value of the order parameter Espi^{Tc)/Tc, when compared with the general bound ( |1.2|) 



confirm the former analyses based on the study of the phase transition: for generic 



5 



We are making here the reasonable assumption that the rate of sphaleron mediated transitions 



in the MSSM is provided by (1.1) and, therefore, that the effect of the new physics in the MSSM is 



entirely encoded in the definition of i^sph- Under these conditions the bound (1.2) still applies 
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values of supersymmetric parameters the phase transition in the MSSM is not strong 
enough to generate the baryon asymmetry of the universe. These negative results have 
motivated a recent search 



J6], |38[ of the region in the parameter space where the phase 
transition has better chances to allow generation of the baryon asymmetry, while being 
in agreement with precision electroweak measurements at LEP [17]. This region has 
been identified as follows: 

• Large values of ttiq guarantee small supersymmetric contributions to the oblique 



radiative corrections 
hancing the value of Rb 



47[| , while small or negative values of can help in en 



48 



Large values of rriA and small values of the stop mixing parameters At and /j, 
favour the strength of the phase transition 



35 



That region was recently analysed in Ref. |^ from the point of view of the strength 
of the phase transition, that was found to be much stronger than in previous analyses. 
Our aim in this section is to study the sphaleron solutions and the sphaleron energy 
in the MSSM, in the above region of the supersymmetric parameter space, in order to 
refine (and eventually confirm) the positive results of Ref. Moreover we will relax 
at some point the condition of large ttia (even though large values of rriA are favoured 
by the strength of the phase transition) because some mechanisms of CP violation 
associated with supersymmetric particles require a sizeable variation of tan/3(Tc) along 
the bubble wall in order to generate the necessary amount of baryon asymmetry at the 
electroweak phase transition pT| . 

The intuitive explanation of why negative values of mfj (keeping the hierarchy 



mh ^ rnfj) favour the phase transition, along with the appearance of color breaking 



4 
minima 

can be found in Ref. 



4^ which can endanger the stability of the standard electroweak minimum, 
]. We refer the reader, for a thorough explanation of this 



phenomenon, to Ref. [^] where all technical details can be found. 

The standard vacuum of the potential ( |2.5|) , with radiative corrections provided by 
(p.8[), has a depth which does not depend on m^, as can be easily checked. Therefore 
the simplest procedure is to compute the depth at the standard vacuum for large values 
of tua'- in this limit the MSSM goes to the SM with a particular value of the Higgs mass, 
determined by the supersymmetric boundary condition of the quartic coupling, and the 
threshold effects due to the stop mixing parameters, at the scale of supersymmetry 
breaking (~ mg). These effects can be encoded in an effective quartic coupling |^ 



and a particular value of the SM Higgs mass, in terms of which the depth of the SM 
minimum is: 

VM = -\^W. (4.1) 



On the other hand, for At = and negative values of ml 



(4.2) 



the potential along the direction f/£ has a minimum at |36 



i m 



u 



(4.3) 
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with a depth 

9s 

Comparison of (|]T) with (4^) allows to determine a critical value of the parameter rnu 



such that the depth of the charge and color breaking minimum ( [4.3|) does not exceed 



that of the standard electroweak minimum ( |4.1j ): 



= [^^yj (4-5) 

For At 7^ we have numerically verified |^ that the condition rnu < "^t/'* also guar- 



antees that the standard electroweak minimum is the global minimum if At ^ 0.8 mg. 
Notice however that such high values of the mixing parameters are uninteresting for 
our analysis since the phase transition will become very weak for those values, as we 
will see. 

We have then fixed mu = m™* and repeated the analysis of Sections 2 and 3. 
We will also fix mq = 500 GeV, to be in good agreement with electroweak precision 
measurements, as we said above, and, for the moment, we will adhere to large values 
of rriA and low mixing: tua = 500 GeV, At = fi = 0, although we will relax the latter 
conditions later on. For these values of the supersymmetric parameters we have plotted 
£^sph at zero temperature as a function of tan/5 in Fig. |10| [solid line] along with the 
comparison, for illustrative purposes, with Ef^ for the value of the Higgs mass given 
by Eq. ( p.29| ) [dashed line], and EsphiTc)/Tc in Fig. |Tl] [solid line], also as a function 



of tan/3. We also plot in Fig. 11 nih [thin dashed line] as a function of tan/5 for the 
corresponding values of supersymmetric parameters. 

Up to here we have used, for our calculations, the critical temperature Tc, where the 
origin gets destabilized and the barrier between the origin and the finite temperature 
minimum disappears. However tunneling with formation of bubbles starts somewhat 
earlier, at a temperature Ti, when the corresponding euclidean action B{Tb) ~ 140 and 
the phase transition can proceed sufficiently fast such that the whole universe can be 
filled with bubbles of the new phase. To compute B{Th) we will use the analytical 
estimate obtained in Refs. [|I^, ^ . First of all we will define an approximated potential 



of the form 

l^app = D{T' - T2)02 -ETct>^ + (4.6) 

where 

(j) = \/2 {hi cos (3 + h2 sin (3) , 

and all the constants have been determined numerically. Tc is the critical temperature 
that we have been using throughout this paper, A = m|/(4t>^), where rrih is the lightest 
Higgs boson mass of the MSSM, 

E = ^ (</>(^c)) 



3 T, 



where {(j){Tc)) is determined numerically, and 



c 
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Secondly, we will use the analytical estimate of Refs. [0, |5^ 



Bin) 



38.8D^/^ /AT\3/2 



f 



2XD /AT 



(4.7) 



where AT = T^ — Tc and the function /, defined in Ref. |5^, equals 1 at zero value of its 
argument. Using (|4.7| ) one can easily determine the value of the transition temperature, 
within the error inherent to the analytical estimate we are using. This degree of accuracy 
is enough for our purposes in this paper. The value of £'sph(Tfe)/Tfe is plotted in Fig. |ll 
[thick dashed line]. From Fig. [T^ we can see that the condition ( |1.2D translates into 
bounds on tan/3 (tan/3 < 3) and on m/j (m/j < 80 GeV). We can also see that one can 
translate the bound ( |1.2| ) at the transition temperature, into a conservative bound at 
the critical temperature T^, 

-E'sph(^c 



T, 



> 50. 



(4.^ 



This bound can be useful in many phenomenological studies of phase transitions. Fi- 
nally, the goodness of the scaling low ( |3.12| ) is shown, for the same values of the super- 



symmetric parameters, in Fig. |T^, where we see it is accurate with an error < 10%. 



To obtain the results of Figs. |I0|-|12| we have been using optimal values of the su- 
persymmetric parameters, from the point of view of the phase transition. However, in 
particular not too large values of rriA might be required, for baryogenesis purposes, as 
we have said above, in order to generate the needed amount of CP violation at the 
bubble walls. In Figs. |13| and we plot i?sph as a function of itla for tan/5 = 2.4, 2.5 
and the same values of the other supersymmetric parameters. In Fig. ^ we compare 
Esph with Ef^ for a Higgs mass mj^ . We see they are very close to each other. In 
genuine two-Higgs situations, when rriA is small, they can depart from each other by 
an amount ~ 1.5%. In Fig. ^ we plot i?sph(^c)/^c as a function of uia- As expected 
we see that the order parameter Esph{Tc)/Tc decreases with increasing uia- Using the 
bound ( [4.8|) puts lower bounds on m^: ttia ~ 100 GeV, for tan/? = 2.4, and > 110 
GeV for tan (3 = 2.5, the corresponding bounds increasing for larger values of tan f3. On 
the other hand, LEP results impose bounds in the (tan f3, rriA) plane. In particular, for 



the set of supersymmetric parameters in Fig. |IJ |51|: niA ~ 120 GeV, for tan/5 = 2.4, 
and rriA ~ 110 GeV, for tan/5 = 2.5, the corresponding bounds decreasing for larger 
values of tan/? |. We can conclude that for tan/5 ~ 2.5, both (|4.8| ) and LEP results 
translate into rriA ~ HO GeV, which is a very safe bound for generation of baryon 
asymmetry in the MSSM [^]. As for the dependence of Esph{Tc)/Tc on the stop mixing 
parameters, we have plotted in Fig. |T3| -Esph (^c)/^c as a function of At, for tan/5 = 1.7 
and the same values of the other supersymmetric parameters. We see that, as expected, 
Esph{Tc)/Tc decreases with increasing At, and the bound (|4.8|) puts an upper bound on 
At as At < 0.4 rnq. 



^We thank P. Janot for sending us some original plots with present LEP bounds. 
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5 Conclusion 



We have constructed the sphaleron solution and computed its energy, at zero and finite 
temperature, in the MSSM for an arbitrary set of supersymmetric parameters. At 
zero temperature we have included the leading one-loop radiative corrections in the 
presence of the sphaleron. At finite temperature we have added the one-loop thermal 
corrections where all daisy diagrams are resummed. At zero temperature we have 
compared -Egph with the corresponding one for the SM with an effective Higgs mass. 
The difference is small 1.5%) even for small values of tjia- We have verified that, at 
finite temperature, the scaling law, where all finite temperature effects can be encoded 
into the temperature dependence of the vacuum expectation value of the Higgs fields, 
is accurate (with an error < 3%), for the cases of weak first order phase transition. 
For the cases of strong first order phase transition (light stop scenario) it is somewhat 
larger (< 10%). Finally our calculation supports the conclusion that for generic values 
of the supersymmetric parameters the MSSM is unable to keep any pre-existing baryon 
asymmetry after the phase transition. However we have confirmed the presence of a 
window for MSSM baryogenesis corresponding to heavy left-handed stops, light (lighter 
than the top) right-handed stops, light CP-even higgses {rrih ^ 80 GeV), not so heavy 
CP-odd higgses (m^i ^ 110 GeV), moderate stop mixing {At ^ 0.4 mg), and small 
values of tan /5 (tan/5 < 3). These results are in qualitative agreement with recent 
non-perturbative calculations performed in the framework of the MSSM |5^, and the 
baryogenesis window will be probed at LEP2. 

Some unconsidered effects should tend to increase the sphaleron energy, and thus to 
improve the previous bounds, while others tend to decrease it and thus to worsen the 
bounds. As for the former, two loop effects at finite temperature have been recently 
proved to strengthen the order of the phase transition in the MSSM [37| and thus they 
are expected to increase the sphaleron energy. As for the latter, first of all, we have 
worked within the approximation g' = in order to consider a spherically symmetric 
ansatz. 0{g''^) effects can be easily accounted for by linearizing the energy functional 
with respect to the U{1) gauge field and neglecting the feed-back on the sphaleron [0, 
or by solving explicitly an axially symmetric ansatz in both cases the SM result 
for the value of the sphaleron energy correction is tiny (~ 1%) as should be in the case 
of the MSSM. A second effect is characteristic of the nature of the MSSM due to the 
large number of scalar fields it contains. We have set all of them to zero in this work []. 
However, in some cases a non-vanishing configuration for these fields might minimize 
the energy functional leading to a modification of the results contained in this paper. 
In any case the variation of E'sph should be negative. In particular, for those fields with 
negligible Yukawa couplings we expect the energy functional to be minimized by the 
zero configuration, since the combined effect of their kinetic terms and the large soft 
breaking masses, which are positive definite, will control the configuration of the energy 
functional. Only in the case of the third-generation squarks in Eq. (|2.27|) , and for the 
case of relatively large values of rnjj and/or At has the energy functional a chance to 
be minimized for non-vanishing field configurations, even in the cases of no charge and 



"^Of course, the equations of motion for all scalars fields of the theory are identically satisfied for 
zero configurations. 
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color breaking at the vacuum. We are at present investigating this issue and the results 
will be published elsewhere |^ . 
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Figure 1: Plots of hi{r)/v, h2{r)/v and /(r) at zero (solid lines) and the critical tem- 
perature T = Tc (dashed lines) for rrit = 175 GeV and the values of supersymmetric 
parameters: tan/3 = 1.5, uia — 100 GeV, mq = 500 GeV, mu — and At — /i — 0. 
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Figure 2: Plot of Esph as a function of tan/? for itlq = 500 GcV, At = fi = and: a) 
ruA = 100 GeV, niu = 0; b) uia = 500 GeV, mu = 0; c) uia = 100 GeV, niu = 400 
GeV; and, d) tua = 500 GeV, mu = 400 GeV. 
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Figure 3: Plot of Egph as a function of rriA for rriQ = 500 GeV, At = ^ = and: a) 
tan/9 = 2, mu = 0; b) tan/5 = 15, m^/ = 0; c) tan/5 = 2, mu = 400 GeV; and, d) 
tan/5 = 15, m?/ = 400 GeV. 
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Figure 4: The same as in Fig. 2 but as a function of the effective Higgs boson mass, 
m'f^. The dashed line is Egph for the Standard Model with a Higgs boson with mass 
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Figure 5: Plot of Egph as a function of At for ttia — rnq — 500 GeV, // = mu — 0, and 
a) tan/3 = 2; b) tan/3 = 15. 
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Figure 6: Plot of Egph as a function of // for uia — rnq — 500 GeV, At — mu — 0, and: 
a) tan/3 = 2; b) tan/3 = 15. 
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Figure 7: Solid [dashed] lines are plots of Esp\i{T) [Espii{0)v{T)/v{0)] as functions of 
T for tuq = 500 GeV, ^ = mu = 0, and: a) tua = 100 GeV, At = 0, tan/3 = 2; b) 
ruA = 100 GeV, A = 0, tan = 15; and, c) tua = 500 GeV, At = 200 GeV, tan /3 = 2. 
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Figure 8: Plot of Esph[Tc]/Tc as a function of tan /3 for the same values of supersymmetric 
parameters as in Fig. 2. 
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Figure 9: Plot of Esph[Tc]/Tc as a function of itia for the same values of supersymmet 
parameters as in Fig. 3. 
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Figure 10: Plot of i^sph at T = as a function of tan/9 for mq = rriA = 500 GeV, 
At = n = Q and mu = m^'\ The dashed hne is the Standard Model value for a Higgs 
mass equal to rrf^ . 
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Figure 11: Plot of £'sph[T]/T, for T = Tc, T5, as function of tan/? for the values of 
supersymmetric parameters of Fig. 10. The dashed line is a plot of the lightest Higgs 
boson mass for the corresponding values of supersymmetric parameters. 
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Figure 12: Solid [dashed] line is £^sph(T) [£;sph(0)i;(r)/'t;(0)] as function of T for tan/3 
2 and mg, m^, ^t, A* ''^u as in Fig. 10. 
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Figure 13: Plots of Esph at T = as a function of tua for tuq = 500 GeV, At = = 
tan/5 = 2.4 (lower solid), 2.5 (upper solid) and mu = m^'*. The dashed lines are 
corresponding Standard Model values for a Higgs mass equal to m^^. 
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Figure 14: Plots of £'sph[Tc]/Tc, as function of for the values of supersymmetric 
parameters of Fig. 13 [upper solid is for tan/? = 2.4 and lower solid for tan/? = 2.5]. 
The dashed lines are plots of the lightest Higgs boson mass for the corresponding values 
of supersymmetric parameters [upper dashed is for tan/3 — 2.5 and lower dashed for 
tan/3 = 2.4]. 
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Figure 15: Plot of Esph[Tc]/Tc, as function of At for tan/3 = 1.7 and niQ, niA, A* 
mu as in Fig. 10. 
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